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We investigate the nonrelativistic reduction of the Bethe-Salpeter amplitude for the deuteron 
electrodisintegration near threshold energies. To this end, two assumptions have been used in the 
calculations: 1) the static approximation and 2) the one iteration approximation. Within these 
assumptions it is possible to recover the nonrelativistic result including a systematic extension to 
OS . relativistic corrections. We find that the so-called pair current term can be constructed from the 

P-wave contribution of the deuteron Bethe-Salpeter amplitude. The form factor that enters into the 
calculation of the pair current is constrained by the manifestly gauge independent matrix elements. 
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I. INTRODUCTION 



In recent years, the discussion of relativistic issues in the reactions involving the deuteron has become more and 
more important. Some well known examples for a clear experimental evidence have been found in the electromagnetic 
disintegration . This results are supported by a recent analysis of polarization observables of the deuteron in a 
relativistic framework Also, the persistent difference between theory and experiment in polarization observables 
of proton deuteron scattering have been argued to be of relativistic origin ||. It is expected that the interpretation 
OS | of new experimental results, e.g., from TJNAF, also calls for a relativistic treatment. 

Meanwhile the number of theoretical investigations devoted to this question and to tackle the relativistic aspects of 
few body systems is increasing. One may basically recognize two lines: on one hand the majority of approaches uses 
a nonrelativistic scheme of calculation and take into account some leading order corrections, such as nonnucleonic 
degrees of freedom like mesonic exchange currents, AA-configurations, A^iV*-states and more. For a review see, e.g., 
Refs. H[|,[i] [| which also contain an large number of references. Note, however, that the problem of consistency arises 
as has been recently pointed out again, e.g., by Ref. jl(| in the context of electrodisintegration. The second line 
^ \ recognizes the importance of covariance JTT|-|T3| . 

We follow the Bethe-Salpeter approach which is covariant by construction and properly denned [p|-|l6| . Based on the 
field theory of particles this approach is also consistent and automatically takes into account all relativistic effects in the 
two nucleon system. Because of the covariance of the Bethe-Salpeter equation the amplitudes (or vertex functions) may 
be written in any reference system, and one may chose a convenient one. However, practical applications are hampered 
by the difficulty of finding solutions of the Bethe-Salpeter equation, even in ladder approximation. Nevertheless, the 
activity on solving the Bethe-Salpeter equation has increased in recent years [|l7|]lq1, and by now parameterizations 



of the Bethe-Salpeter amplitudes in Euclidean space are available |19[|. Also, useful results have been achieved by a 
reduction of the Bethe-Salpeter equation to a three dimensional one. 

Here, we examine the nonrelativistic reduction of the Bethe-Salpeter approach to compare with the results achieved 
by nonrelativistic calculations. Since the problem has not been solved so far for the general case, we restrict ourself 
to a special process that is the disintegration of the deuteron. To this end we consider the threshold region only that 
exhibits already the main features to be discussed. Also, this region is dominated by a single transition amplitude to 
the 1 S'o final state (see, e.g., Ref. ||(|) thus leading to some technical but not substantial simplifications. The extension 
to other partial amplitudes is straightforward but tedious, however beyond the scope of our present approach. 

Indeed, the deuteron disintegration close to threshold energies has been a key reaction to establish subnucleonic 
degrees of freedom in nuclei |2M and it attracts continuous attention from theoretical and experimental sides. It has 
been an excellent paradigm to examine nonnucleonic degrees of freedom and relativistic effects. It is well known that 
the nonrelativistic impulse approximation fails to describe the double differential cross section at momentum transfer 
squared — q 2 > 9 fm~ 2 . The experimental data do not indicate the deep minimum present in the calculations [p2[ . 
To fill this minimum mesonic exchange currents have been introduced. The contributions of 7r-, p-currents, and AA- 
configurations allow one to get a satisfactory agreement with the data |22j ]. The question of a consistent inclusion 
of all relativistic corrections (at least for the pion sector) has recently been addressed by Ref. (Tfl] and supports the 
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above statement. Nevertheless, some conceptional problems of the theory still remain open questions. Among them 
is the problem of gauge invariance and the question of the nucleon form factor to be chosen in the exchange currents 
contributions |24|^2q|. 

By now some relativistic methods identify corrections like meson exchange currents in the relativistic impulse 
approximation. It has been shown that in the framework of the light-cone approach the so-called extra components 
fs of the deuteron wave function, and g 2 of the 1 So-wave function introduced in Ref. fl27|] give an expression in 
the nonrelativistic limit, that analytically equals to the contribution of the pair current p7|. In this context the 
one iteration assumption that will be explained below for solving the dynamical equation is substantial. Another 
important result in this context is the calculation of the static electromagnetic properties within the Bethe-Salpeter 
approach. It follows that the contribution of the P-states to the magnetic moment of the deuteron is numerically 
close to the contribution of mesonic currents and agrees in sign p8| . It has already been noticed earlier in the context 
of covariant reductions of the Bethe-Salpeter equation to three dimensional ones that negative energy components 
in the wave function are responsible for pair current type contributions, see, e.g., Ref. [ p9| . Here we show that the 
suggested reduction procedure leads to the same analytical structure as the nonrelativistic pair current correction. 

The paper is organized as follows. Section [n] contains the general formulas for the deuteron disintegration amplitude 
in different rep resentations. The static approximation is also discussed. The one iteration approximation is formulated 
in Section |BJ where we give reasons for its validity in simple models. The derivation of the main correction is done 
in the Section |iy| We present and discuss our results in Section [v|, summarize and conclude in Section VI . 



II. NONRELATIVISTIC REDUCTION OF THE ELECTROMAGNETIC CURRENT MATRIX ELEMENT 

Because of Lorentz covariance, transformation properties under parity and time reversal as well as current conser- 
vation, the general structure of the transition matrix element for the 1 + — > + transition current near threshold can 
be written as |30| 

{npCS )\j^\DM) = ienafr q V(s, q 2 ), (1) 

where is the deuteron polarization four vector with projection A4, K is the deuteron four momentum, q = P — K 
is the four momentum transfer, s = P 2 = (K + q) 2 w 4m 2 is the squared total momentum of the np-pair, see Fig. [I]. 
The scalar function V(s,q 2 ) may be splitted into isovector Dirac (F^) and Pauli (F^) contributions 

V(S, q 2 ) = V ± {S, q 2 )F[ V) (q 2 ) + V 2 (s, q 2 )F^ V \q 2 ). (2) 

The final state of the np-pair is supposed to be in a 1 So-state. This is supported by nonrelativistic calculations |p3| . 
Eq. (Q) is valid in the any reference system and the functions depend on the scalars s and q 2 only. Current 
conservation is fulfilled because the tensor e A1Q/ 3 7 that appears in Eq. dTl) is antisymmetric. 

To calculate the functions Vi,2 from the underlying dynamics in the Bethe-Salpeter approach, two representations 
of the Bethe-Salpeter amplitude may be used in concrete calculations: 1) the covariant representation and 2) the 
partial-wave representation. Using the covariant representation for the initial and final states allows us to directly 
read off the functions Vip(s,q 2 ) from the proper current matrix elements. This has been done explicitly in Ref. [ p0[ . 
However, to study the relation to the nonrelativistic expressions it is more suitable to give the formulas for Vi^{ s >?) 
in terms of a partial-wave representation of the vertex functions. Hence, the expressions for Vi,2 can be written as 

V h2 (s,q 2 )=i J d 4 k^2 G a (p ,\p\) O«0 9p(ko,\k\), (3) 

where the partial vertex functions G\-4,{j>q, \p\) represent the states 1 5 , f+ , 1 S ~~, 3 Pq, 3 Pq of the np-pair, and 
gi- 8 (k , \k\) denote the states 3 S' 1 ++ , 3 S7~, 3 D^ + , 3 Dj_~, 3 Pi, 3 P?, 1 P°, 1 P° of the deuteron. We use spectroscopic 
notation and (±, e/o) as p-spin quantum numbers |31[ . In general, the lengthy functions O a p depend on Lorentz 
scalars and the explicit expressions are omitted here. We will specify them below after having introduced appropriate 
approximations. The expressions for po and |p| are given in a formally covariant way, viz. 

n _ I* + </.* + 0/2) , w _ + + _ (k + ?/2) ,) (4) 
where ( , ) denotes the Lorentz scalar product. 
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We note that through the use of Bethe-Salpeter vertex functions the denominators of O a /3 appearing in Eq. (j^) 

contain products of (±M/2 ± ko ± E k ) and (±y/s/2 ± po ± E^), where _E p = \J p 2 + m 2 , etc. that stem from the 
nucleon propagators (see Fig. |l|). We evaluate the integrals in the laboratory system (deuteron at rest). At threshold, 
because of the small deuteron binding energy it is possible to utilize the static approximation f32|| that preserves the 
analytical structure of Eq. (m through the following equations 



/>« = /.■„. \P = + - (k 2 + ^ + \k\\Q\x) . (5) 



Here and in the following we use (a = a/|a|) 



Q =\Q\ = {{M 2 + s - q 2 ) 2 - 4Af 2 s) 1 / 2 /(2M), x = k-Q, 

q= (u,Q), oj= (s- M 2 ~q 2 )/(2M). (6) 

We illustrate the approximation by looking closer to the expressions involving the (++) states. The full denominator 
of the integrand in Eq. ([|) leads to a complicated pole structure and reads explicitly 

^2 i „2 rfQ 2 /i , u \2 , r 2\ 1/2 , .\ ( M 



^4-^0+2) Z +P -Vs[^-(k + -y + E^j -iej^ Y + k -E k + iej. (7) 
Using Eq. (0), then leads to a simple pole expression for the integrand involving the denominator 



Js \ ( J~s \ (M 



ko- Ep + ie) ^- -ko- E p + ie) — + fc - E k +ie) . (8) 



Thus, static approximation means in particular that u) = (no retardation) and the Lorentz boost transformation of 
the np-pair vertex functions is neglected. Going beyond the static approximation can be achieved by expanding the 
full expression in terms oicu/M and Q 2 /s that lead to additive corrections. 

The integration on kg can now be performed by choosing a proper integration contour and specifying the corre- 
sponding poles, e.g., closing the upper half plane leads to poles for kg at po — y/s/2 — E p . The vertex functions are 
then evaluated at (po,k). Since in the reaction under consideration s rs 4rn 2 ss M 2 , we expand the vertex functions 
near ko = M/2 — E k for g\ and G\ respectively that allows us to derive analytical expressions in the one iteration 
approximation as will be shown in the next section. The analogous procedure holds for the other partial vertex 
functions. With this choice one the nucleons in the deuteron is taken on-shell. 

Alternatively it would have been possible to expand the functions, e.g., near kg = 0. However, the above advocated 
approach is more suitable, since it leads to an equation that allows for an analytical of /co- 
in Eq. (^) we now perform the kg integration as explained above. The angular integration is simplified by taking 
Q along the z axis and replace x given in Eq. (|^) via 

k /2 -k 2 -Q 2 /4 , 2k' „, 

where k = \k\. Finally, to compare Eq. (|J) with the nonrelativistic limit we introduce a /c/m-expansion excluding 
order 0(k 2 /m 2 ) (i.e. E k = m + Oik 2 /m 2 ) — E p ). The resulting structure functions V\.2(s,q 2 ) are then given by 

+ oo k+Q/2 . 

Vl {a,f) = ^-± [ kdk I k'dk'l ^ k ' ] ( gl(fc °' fc) 4^ P 2 (x) 
u ' mMQ J J \Vs- 1E k , +ie\M- 2E k \f2 M — 2E k 2y ' 



\k-Q/2 
V3 Gi 



g 5 (k 0l k) -g 6 (k ,k) 



-G 3 (k' ,k') + G 4 (kik') 



2 ^/s-2E k ,+ie 

V2 t 

4 

.9i(fc ,fc) Q + 2kx 1 g 3 (fco,fc) QP 2 (x) + 2kx " 

M - 2E k Qk' 71 M-2E k Qk> , 



(10) 
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+ oo fc+Q/2 . _ 

2 ^ ,q ' ~ mM Q J J kdh \ VS-2E k ,+ie\M-2E k M - 2E k 

\k-Q/2\ K 

, V3 Q 2 Gi(^,fc') 



4 4™^-2E, +K ^ fe ^ 9 *' fc) jQ 
3^2 Q 2 



gi(fco,fc) Q + 2fca; _ J_ g 3 (fc ,fc) Qf 2 (i) + 2fcx 
M - 2S fe Qfc' ~ ^2 M - 2E k Qk' 



(11) 



where k' = p = t/s/2 — E k > , and fc = M/2 — E k , and Pa^) = (3a; 2 — l)/2 is the Legendre polynomial. The functions 
92, 54) 97, 58: Gi disappear in the above expressions after fco integration and because of the k/m expansion. Note, 
that within this approximation scheme we are left with the (++) to (++) and (++) to (e, o) transitions only. All 

other matrix elements, such as ( ) to all and (e, o) to (e, o) cancel, respectively. 

We now examine the expressions for Vi^{s,q 2 ) given in Eqs. more closely. To recover the nonrelativistic 

result, we neglect the vertex functions g$, ge, G3, G4 that correspond to the negative p spin components (i.e. do 
not exist in the nonrelativistic scheme). If we replace the functions g\ and g% by the nonrelativistic S and D wave 
functions and G\ by the 1 Sq continuum wave function in the following way 



9i(kp, k) 
M — 2E k 

53(fco,fc) 
M - 2E k 



-a±u(k), a\ = 47iV2M 



-cti w(k), 



GiCK,k') ,,,, 11 . 

a 2 u (k), a 2 ^—i= — , (12) 



Vs - 2^ fe / + is V47T 27T 

and insert these for the respective vertex functions into Eqs. (JlOjll) we obtain 

+ 00 fe+Q/2 

v(0){s > q2) = ^rh G ** ){q2) I kdk I k ' dk ' u o( k '){ u ( k )-^ w ( k ) p ^))- ( 13 ) 

|fc-Q/2| 

where we have introduced the magnetic isovector form factor — + F 2 V ^ . In co-ordinate space the respective 
integral is achieved using the following transformations, for the deuteron states (w = u, w 2 = w, I = 0, 2) 

+ oc +00 

wi(k) = J rdrwi(r)j t (kr), ~~~ = ~ / k 2 dkwi(k) j e (kr), (14) 


where ji(x) is spherical Bessel function, and for the scattering state 

+oc +00 

Uo(k) = — J rdru (r)j Q (kr), °^ - = / k 2 dku (k) j (kr). 



The resulting expression is 

v[0){s ^ )= a^ G {v) {q2) J dru (r)(u(r)j (Qr/2)--Lw(r)j 2 (Qr/2) S ), (15) 



This result reflects the so-called nonrelativistic impulse appro ximation and represents the lowest order nonrelativistic 
expansion of the transition form factors given in Eqs. ( Eo| , ). 
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III. ONE ITERATION APPROXIMATION 



A. Deuteron channel 



The Bethe-Salpeter equation is commonly solved by iterations. After angular decomposition we are left with an 
integral equation for the radial parts of the vertex function g a {po,p)- They are connected to the partial amplitudes 
if>0{ko, k) which are used here for simplicity on the right hand side of the following equation by amputation jjlj, e.g., 
for the (++) components the relation is given by 

^ = (M/2 + k -E^(M/2-k -E k y a = 1 ' 3 ' (16) 

The integral equation read 

+00 +00 

9a(po,p)= / dfc ° / ^^^dkJ2vlf\po, P ;k ,k)Mko,k), (17) 

M -00 p P 

where M denotes the type of exchanged meson, gMNN its coupling constant and V^gP is the transition matrix element 
between the states a and (3 given below. To get fast convergence to the solution of this equation one needs a good 
educated guess for the initial vertex function g a {Po,p) respectively ^(^cb k). This guess may be taken from the 
solution of the equivalent nonrelativistic problem. After several iterations one usually gets the exact solution. 
The initial vertex functions may be chosen along the lines given in Ref. |^], generalizing Eq. ( [l2| ) 

,, (h ,x_ -<*i (M-2E k )u(k) 

Wo, K) - (M/2 + kQ _ Ek+ ie){M/2 _ h) _ Ek+ i6 ) > W 
h(h h s_ -"1 (M-2E k )w(k) 

Wo, K) - (M/2 + kQ _ Ek+ i£)(M/2 - ko -E k + ic) ' 1 ' 

All other wave functions are assumed to vanish in this order. Note, that the numerators of Eqs.( |T^Jl9| ) are the 
nonrelativistic vertex functions. 

The ansatz given in Eqs. (|l^) and ([l9|) has been numerically verified for the deuteron consisting of two spinor 
nucleons that interact via the exchange of scalar mesons with mass \x = 0.2 m. This interaction leads to a bound 
state for both, nonrelativistic and relativistic cases. In the nonrelativistic case and use of the Schrodinger equation 
the kernel is of the form 

v{p - k) = 9jm — _L — (20) 

v ' Air (p - k) 2 + n 2 y ' 

where the coupling constant gNR, — 8.37 leads to a binding energy of e — 2.23 MeV. The resulting wave function of 
the deuteron may be parameterized in the following form, 

_ / 1 0.40 0.60 \ 

U[P) ~ \WT^ ~ p 2 + (A + fi) 2 ~ p 2 + (A + 2/i) 2 ) ' { } 

with A = y/me and a normalization constant N. Note, that for the scalar exchange model the I?-wave function is 
zero, w(p) = 0. These functions are then substituted as an educated guess into the system of equation ( |l7j ) using the 
assumption of Eqs. ( ^8p9| ). Eq. (|l7|) is then solved in Euclidean space (p — > ip±) using the ansatz amplitude Eq. ( |2l| ) 
and the full kernel V(pq,p; ko, k). 

a 2 1 

V(p ,p; ho, k) = |- _ _^2 1 ® 1 ( 22 ) 



To achieve the same deuteron binding energy as before, now using the Bethe-Salpeter equation (|17|) the coupling 
constant g of the scalar interaction needs to be readjusted to get the same deuteron binding energy, viz. g = 17.00. 

In this framework we compare the full vertex function of the exact solution of Eq. (JT^) to the result achieved 
by one iteration only but using the nonrelativistic solution as an educated ansatz. This will be called one iteration 
approximation in the following. As a demonstration for the quality of the approximation we show two vertex functions, 
for the 3 Si + - and 3 P°-channels in Fig. ||. The upper two surfaces displayed correspond to the exact result and the 
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one iteration approximation of the 3 5^ -component. The difference between the two is rather small. In fact, in 
the relevant region p = 0. . .0.5 GeV and p^ = —0.1 . . .0.1 GeV, where the amplitude has its major contribution 
to the normalization and to physical processes, the difference does not exceed 5%. The exact and the one iteration 
approximation for the 3 P 1 °-channel that is relevant for the subsequent discussion of the electromagnetic properties 
are displayed as the lower two surfaces. We conclude that for the scalar model the shape of the exact solutions are 
already reproduced by one iteration. From that we argue that the one iteration approximation is a reasonable first 
step that may be utilized to solve the more complicated problem with realistic potentials. 

We now turn back to the discussion of the realistic deuteron implementing the full kernel with all spin dependences. 
By using Eqs. ( |i~7j ) along with Eqs. ( PlP^t ) the expressions for the 3 Pf ( 3 P°)-channel, i.e. for the functions g${go) are 
now deduced analytically. If we consider the 7r-exchange kernel only and use the above given notation for the states 
(1...8) then 

V 51 (p ,p;k ,k) = V 53 (p ,p;k ,k) = 0, (23) 
V 61 (p ,p;k ,k) = -j= (^p(E k -m)Q 2 (z)+mkQi(z) - (E k + 2m)Q {z^J , (24) 

v/2 / 1 2 \ 

V 63 (po,p;ko,k) = — i-p(2E k +m)Q 2 (z) ~mkQ 1 (z) - -p(E k - m)Q (z)J , (25) 

with Legendre functions of the second kind Qe(z), and z — (p 2 + k 2 + fj, 2 — (ko — po) 2 )/2pk. The ko integration can 
be performed yielding (e.g. for the pionic exchange kernel) 

9aiP °' p) = £^^ir I -^{V^{p,k)u{k) + V^{p,k)w{k)), a = 5,6. (26) 



The functions V^j 3 \(p, k) are obtained from V^j 3 \(j>o-,P', ko, k) by treating po and fco as explained in Sect. ||. The 
explicit expressions for V (5, 6 <— 1, 3) are given by: 

V 5 i( P ,k) = V 53 (p,k) = 0, (27) 
2m 

V 61 (p,k) = -j=(kQ 1 (z)-pQ (z)), (28) 
\/2m 

V 63 (p,k) = ^-=-(pQ 2 (z)-kQ 1 (z)), (29) 

with z = (p 2 + k 2 + /j, 2 )/2pk. To simplify the notation the isospin factors have been omitted in V and Eq. (p6|). The 
extension to other mesons of a one boson exchange kernel is straight forward. 

Now using Eqs. (|2^,|2^), the transformation Eq. (|l4| ) and the integral representations of the Legendre functions 

+ OG 

Qe(z) — 2pk J e~ tJ,wT r dr je(kr) ji(pr) we obtain following results for the vertex functions: 
o 

g 5 (k ,k) = 0, (30) 

-j-OG 

g 6 (h,k) = -(-3)^-^. [ dr e —{l + ^r)(u{r) + \w{r)]ji{kr), (31) 



V3m Air J r V V2 

o 

where the isospin factor (—3) has now been shown explicitly. This is the analytical structure given after the first 
iteration. 

B. 1 So-channel 

The inhomogeneous Bethe-Salpeter equation for the amplitudes in the ^o-channel reads 

+ OG +00 







x S a ff(p ,p; s)V^ '(po,p;k , fc)0 7 (fc o , fc). (32) 
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Here (j)^ denotes the plane-wave function, 



"(po,p) = 4= S al S(p ) \ S(p-p*), (33) 

V47T V 



where p* is the on-energy-shell momentum given by p* = \p* \ = \J s/A — m 2 . 

For the subsequent discussion it is more convenient to split the system of equations into the following form: 



,-oo +OC 

-,2 



4£ fc £ p P 

y^ V i0 (Po,P;k ,k) , , (34) 

^ p (Vs/2+po - E p + ie)(Vs/2 -p - £ p + ie) 



M 
X 



X) v i?(Po.P;*d,fc)^(*b,*). o , ^l= 1 S+ + . (35) 

Since the full vertex functions are related to the full amplitudes by amputation, the partial vertex functions G a {po,p) 
are connected to the partial amplitudes <fi a (j>o,p) via simple relations As in the case of the deuteron channel we 
consider the one iteration approximation. To do so we chose a similar expression for the initial amplitudes as before 
to connect to the nonrelativistic solution, 

, (h .A -a 2 {Js-2E k ,)u Q {k) 

M k , ft ) - (y _ /2 + ^ _ ^ + 2£)(y _ /2 _ ^ ^ + . e) . (3b) 
Here uo(k) is the nonrelativistic continuum wave function in 1 5'o channel given by 

no(fc) = F ^~p)+ p * 2 _ fc2+ , £ > (37) 

and t(k,p*; E p *) is nonrelativistic half-off-shcll i-matrix for the 1 Sq channel normalized through the condition 

t (E p , ) = tip*, p*;Ep.) = _ 1J_ sin <J e w ° , (38) 

ix rap* 

where (5o is the phase shift, and Sp* = p* 2 /m. 

Analogously to the deuteron case, we finally arrive at the first order corrections to the amplitudes in the channel, 

G 3 (fco,fc') = -(+l)^4 m / dr e —-{l + ^r)uo{r)h{k'r), (39) 
irm 47r J r 
o 

G 4 (k ,k')=0, (40) 

where (+1) is the isospin factor. 

We have shown in this section that choosing proper zero approximation wave function (i.e. the nonrelativistic ones) 
after one interaction additional partial amplitudes arise through the Bethe-Salpeter equation. They are connected 
to the interaction kernel and in electromagnetic processes give rise to the so called pair current correction as will be 
shown in the next section. 



IV. THE STRUCTURE OF THE ELECTROMAGNETIC CURRENT MATRIX ELEMENT 

We are now in the position to turn to the first order corrections to V^°\s, q 2 ) given in Eq. (|l5|). To this end 
we expand the expressions given in Eqs. (g), (|lfj) and (|TTJ) into a power series of g\ NN j\-n to extract the pionic 
contribution only. Also we consider the P-states contribution only, i.e. components with one negative p spin. The 
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resulting expression for the transition form factor V(s, q 2 ) will be denoted by V^(s, q 2 ). Substituting Eqs. (po| ) and 
( |3l"l ) as well as Eqs. ( j39[) and (|io| ) into Eqs. (|o|) and ( |ll| ) and using the replacements of Eq. (|l2|) then Eq. (g) reads 

+ 00 +00 fc+Q/2 

y«( S) ^ ) = ^I^%iv f dr ^l {1 + ^ r) f kdk f y dk , 



2m 2 MQ 4tt 

\k-Q/2\ 

x|(-3) (VfV)^^fV)) uo (*o(«(r) + -^«;(r))j 1 (fcr)| 
+(+1) (^V)-^^ 2 )) i UoWjl (fcV) 



, Q + 2fcz 1 QlMx) + 2kx\ 
X \ u ( k )—Fnj 7^ w ( k ) If- ( 41 ) 



Qk> ^2 Qk' 

The fc'-integration can be solved analytically. After doing so, we obtain 

+00 

yM {oa 2 ) _ a ^ 91nn f dr ^l a + ur ) 
V ^ ' 2m 2 MQ An J r 2 ^ + ^ r > 



x| 3 (*fV) - \^2 F 2 V \l 2 )) Mr) (u(r) + ^f(r)) Ji(Qr/2) 
ifV) ~ l^2 F 2 V) (l 2 )) Mr) («(r) + ^u(r)) ji(Qr/2)| 

= ^^ H ^ 2 )% L / ^^(l + ^uoW^W + ^W^^Qr^). (42) 


Here we have introduced the function 

^ 2 )=^V)-^FfV) (43) 

where Q is given in Eq. (^|). 

This first order contribution in g 2 NN supplements the lowest order relativistic expansion as given in Eq. (|l^) of 
Sect. ||. We then arrive at the following expression for the transition form factor: 

+00 



V(s,q 2 ) = ^^lG^(q 2 ) J dru (r) (u(r) j (Qr/2) - ±w(r) j 2 (Qr/2) 



+00 

On 2 f P ~fJ-7,r / I 



+i/(<?) ^ ™ J d r-^(l + ^r)u (r)\u(r) + -^w(r))h(Qr/2)\ (44) 
' ' 

This is our main result. Starting from a fully covariant theory this result has been achieved using the lowest order 
relativistic contribution within the one iteration approximation. The relation between the nonrelativistic wave func- 
tions and the relativistic amplitudes given in Eqs. (|l9|), and ( p6| ) allows us to give the final expressions in terms 
of nonrelativistic wave functions. The lowest order in the irNN coupling constant g 2 NN leads to an additional contri- 
bution from iterating the P-wave channel once. Comparing this result to the one achieved within the nonrelativistic 
scheme that introduces mesonic exchange currents, we find find that the first term coincides analytically with the 
nonrelativistic impulse approximation contribution and the second one with the 7r-pair current contribution. 



Nonrelativistic impulse approximation and pair current contribution 

To further illustrate the equivalence we connect our results to the nonrelativistic formalism describing the break-up 
process. More details are given, e.g in Ref. |Q. The differential cross section for the D — > 1 Sq - transition has the 
form (nonrelativistic case) 
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d l G 16 2 k 12 p*m 



-a 



dn<kj - 3 ~ Q2 ^ S in 2 -(( fce+ ^-2 fc ^cos 2 -)|(%||Tr E ||D)^ (45) 

where q = (u>,Q) is the momentum transfer, t — —q 2 . The momentum p* is related to the relative energy E p * 
of the np system as given before E p * — p* 2 /m and the relation between kinematical quantities is given by Q = 

A/((2m + E p ,) 2 - M 2 + t) 2 /4M 2 + t and k' e = {-lu + yju 2 + f/shr 2 0/2)/2, w = P e - ££. 

In the general case, the current matrix element is a sum of nonrelativistic impulse approximation contribution, 
meson exchange currents contributions and retardation currents contributions. Here we are interested in 7r-meson 
pair current part only and thus we may write for { 1 Sq \\ T^ as || D) the following expression: 

CS || lf ag || D) = ( l S II || D) + CS || || D) (46) 

where we introduced Pf^ 8 , which reflects the impulse approximation operator and Tj^f , which is the n- meson pair 
(contact) operator. Expressions for the matrix elements are given in Ref. fl34| , see Eq. (11), and may be compared 
with Eq. Differences will be discussed below. 



V. RESULTS AND DISCUSSION 



We have shown in the previous section that the nonrelativistic reduction of the Bethe-Salpeter approach utilizing the 
one iteration approximation leads to results that exhibit the same analytical structure as the nonrelativistic result plus 
pair current corrections. Some details differ as will be discussed below. One may now use the "exact" nonrelativistic 
wave functions to calculate the different contributions to the cross sections. This has been done here for an illustration. 

The formula for the cross section is given in Eq. ([l5]). The dominant Ml transition matrix element (i.e. D — > 
np^So)) of the multipole decomposition given there directly corresponds to the nonrelativistic reduction of Eq. (|l|) 
given above. The contribution of the nonrelativistic impulse approximation given, e.g. in Ref. |m| coincides with the 
formula given in Eq. ( ]l5| ) if ct\ and a.2 are chosen as in Eq. (12). 

The analytical structure of the nonrelativistic pair current contribution equals to that of the P-state contribution 
derived from the Bethe-Salpeter approach given in Eq. (H3). We note that the nonrelativistic pionic pair current 



contribution given in Ref. ||34j depends on the nucleon form factor Fy{q 2 ). Sometimes also the electric nucleon 

Sachs form factor G ( p(q 2 ) = F 1 (V) {q 2 ) + j^F 2 [V) {q 2 ) has been used (see Ref. @). The Bethe-Salpeter approach 
yields a different dependence on the nucleon form factor that is also consistent with current conservation and given 
by the function H(q 2 ) defined in Eq. ([43|). As an illustration of the different behavior we display the form factors 
F [ l V \q 2 ) (solid line), G ( p (q 2 ) (dotted line), and the function H{q 2 ) (dashed line) in Fig. | As a parameterization 
for the nucleon form factors we use the one given in Ref. |35|]. The form factor F\{q 2 ) is larger than the form factor 
GE{q 2 ) and the function H{q 2 ) and therefore the respective pionic pair current is expected to be larger than the one 
using the other form factors. Since the function H(q 2 ) is in the middle of the other ones the respective contribution of 

the P-states in Bethe-Salpeter approach will be different from the nonrelativistic calculations normally using Fy {q 2 ) 
orGgV); 

To investigate the influence of the nucleon form factors more closely we calculated the impulse approximation and 
pionic pair current contributions to the differential cross section. The calculation is performed with the Paris NN 
potential |§9| at E np = 1.5 MeV and = 155°. 

It is well known that some uncertainty is related to the hadron form factors. Without dwelling too much on 
that point we would like to compare three different form factors to see how this uncertainty influences the result. 
Introducing the hadron form factors change the expressions for the pair current contribution p8[ . Three sets of hadron 
form factors (for it N N- vertex) have been utilized and are shown in Figs. ^ a-c), respectively: a monopole vertex |36| ] 
with a cut-off mass of A = 1.25 GeV (set a) and A = 0.85 GeV (set b). In addition a vertex inspired by a QCD 
analysis has been used with two parameters given chosen to be Ai = 0.99 GeV and A 2 = 2.58 GeV |57]] (set c). 

It is seen from Fig. ||that there is a strong dependence of the differential cross section on the nucleon electromagnetic 
form factors. It is also seen that the minimum at t = 12 fm~ 2 in the impulse approximation contribution is shifted 
by pionic pair current to the region t > 16 fm~ 2 using Ge((7 2 ), to t > 18 fm~ 2 using the function H(q 2 ) and to t > 22 

fm~ 2 using Fy V \q 2 ) as a form factor. The largest shift in the cross section is achieved by using F[ V \q 2 ) in the 
calculations. The size of the shift and the behavior of the cross section considerably depends on the set of parameters 
as well as the type of the irNN- vertex used. This is illustrated in Fig.||d) for calculations using Fy(q 2 ) only but for 
different parameterizations of the strong vertex. 
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VI. SUMMARY AND CONCLUSION 



We investigate the electromagnetic current in the framework of the Bethe-Salpeter approach. The question ad- 
dressed is the nonrelativistic reduction of the current transition amplitude. As a specific application we consider the 
deuteron disintegration near threshold energies in some detail. We present a method to reduce the amplitude to the 
nonrelativistic one. We summarize the basic steps of this method and the main approximations utilized. 

To arrive at our result, i.e. to compare the relativistic expressions to the pair current that appears in the non- 
relativistic scheme, two main assumptions have been introduced: 1) the static approximation and 2) one iteration 
approximation. In the static approximation boost corrections are neglected. We emphasize that the latter contribu- 
tions are not negligible and for the elastic form factors corrections can achieve 20 %. A discussion of this contributions 
in a noncovariant but consistent relativistic framework may be found in Ref. |^0| . Although important, a discussion 
of all these contributions in not the main issue of our present paper. In fact, we do not expect any substantial change 
in our present conclusions, since all corrections are of additive nature and will be taken into account as the scheme 
develops. 

The one iteration approximation is to some extent more subtle. It plays an essential role in connecting the Bethe- 
Salpeter amplitudes and the usual nonrelativistic wave functions. This approximation allows us to express the final 
result in terms of the usual S- and I?- wave functions of the deuteron and the S- wave function of np-pair. Although 
numerically verified for deuteron-like models, the one iteration approximation needs further investigation. We argue 
that the validity is due to the small binding energy for the deuteron. From inspection of the next iterations we find 
additive corrections only. 

Therefore these assumptions can be loosen up and additional corrections can be taken into account. This way we 
have a possibility to calculate corrections in a more consistent way. This problem is present in both, electromagnetic 
and strong processes (for pd — > pd and pd — > pX reactions, for instance). 

Although the calculation of these corrections is important to arrive at definite conclusions about the covariant 
approach as a whole, our main issue here is to investigate the correspondence between the covariant approach and 
the mesonic exchange currents of the nonrelativistic approaches. We have shown the analytic correspondence of these 
different approaches. 

The pair correction plays an essential role in explaining the experimental data. An interesting result in this context 
is connected to the form factor entering in the pair current term. The overall form factor that is consistent with 
gauge independence and the relativistic (covariant) expansion neither coincides with Ge nor with F\ but requires an 
"average" between them. 
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FIG. 1. Impulse approximation of the deuteron disintegration process. 
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2 rr -2, 

t--q [fin [ 

FIG. 3. Nucleon electromagnetic form factors for the nonrelativistic pair current contribution as discussed in the text. Solid 
line F[ v \q 2 ), dotted line g£°(? 2 ), and dashed line H{q 2 ) as given through the Bethe-Salpeter approach. 
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FIG. 4. Impulse approximation (ia, solid lines) and 7r-meson pair-current contributions (pc) of the differential cross section, 
a-c) different sets of hadronic form factors as explained in the text. Dashed line: ia+pc with F^ v \q 2 ), dotted line: ia+pc with 
H(q 2 ), long dashed line: ia+pc with G^\q 2 ). d) Calculations using F^ v \q 2 ) only but different hadronic form factor sets as 
given in the text. Dashed line: ia+pc with set a, dotted line: ia+pc with set b, and long dashed line: ia+pc with set c. 
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